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Abstract 

Investigations  for  the  melt  growth  gallium- an timonide  (GaSb)  semiconductor  crystals  are  underway  at  the  US  Air  Force 
Research  Laboratory  at  Hanscom  Air  Force  Base  by  the  vertical  gradient  freeze  (VGF)  method  utilizing  a  submerged  heater.  Elec¬ 
tromagnetic  stirring  can  be  induced  in  the  gallium-an timonide  melt  just  above  the  crystal  growth  interface  by  applying  a  weak  radial 
electric  current  in  the  melt  together  with  a  weak  axial  magnetic  field.  A  force  convection  in  the  melt  can  be  induced  by  rotating  the 
heater.  This  paper  uses  a  Chebyshev  spectral  collocation  method  and  investigates  the  effects  of  heater  rotation  and  electromagnetic 
stirring  on  the  melt  motion  during  the  VGF  process  using  a  submerged  heater.  Clockwise  and  counterclockwise  rotations  of  the 
heater  augment  and  opposes  the  flow  induced  by  the  electromagnetic  stirring,  respectively.  The  radial  velocity  along  the  crystal-melt 
interface  is  proportional  to  both  the  strength  of  the  electromagnetic  stirring  and  the  clockwise  rotation  of  the  heater.  A  relatively  low 
rotation  rate  in  the  clockwise  direction  significantly  increases  this  radial  velocity. 

©  2005  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Bulk  gallium-antimonide  (GaSb)  semiconductor  crys¬ 
tals  with  high  optical  transmission  are  extremely  im¬ 
portant  for  space-based  imaging  applications.  These 
crystals  can  be  grown  from  the  liquid  phase  (melt)  and 
the  melt  is  typically  doped  with  an  element  in  order  to 
give  the  crystal  certain  optical  properties.  The  growth 
of  these  crystals  from  the  melt  is  critically  important  be¬ 
cause  viable  substrates  which  are  compositionally  uni¬ 
form  both  within  a  wafer  and  from  wafer  to  wafer  are 
needed. 

Bulk  gallium-antimonide  crystals  can  be  grown  from 
the  melt  by  the  liquid-encapsulated  Czochralski  (LEC) 
process  or  by  the  vertical  gradient  freeze  (VGF)  process. 
The  LEC  process  is  a  top-seeded  growth  method  that 
has  distinct  advantages  because  the  crystal  is  viewed  di- 
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rectly  during  growth  and  because  there  is  no  container 
in  contact  with  the  crystal.  However,  short-range  com¬ 
positional  uniformity  can  be  problematic  due  to  buoy¬ 
ancy-driven  convection  in  the  melt,  which  randomly 
changes  the  composition  of  the  melt  in  the  boundary 
layer  at  the  growth  interface  (Kuniholm  and  Ma, 
2003).  In  addition,  long-range  crystal  composition  gen¬ 
erally  changes  due  to  a  changing  melt  composition  for 
dopants  with  segregation  coefficients  not  equal  to  unity 
(Morton  et  al.,  2002).  Relatively  good  long-range  com¬ 
positional  uniformity  was  achieved  for  a  non-stoichiom- 
etric  Sb-rich  melt  by  the  LEC  process  by  Ohmori  et  al. 
(1982).  However,  only  25%  of  the  melt  volume  was 
solidified  because  the  continual  increase  of  antimony 
in  the  melt  as  crystal  growth  progressed  eventually 
caused  the  crystal-melt  interface  to  breakdown  due  to 
constitutional  supercooling.  In  the  early-90s,  Ostrogor- 
sky  (1990)  introduced  a  modification  of  the  bottom- 
seeded  VGF  process  in  which  a  submerged  heater 
separates  the  melt  into  two  zones,  namely,  a  lower  melt 
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and  an  upper  melt.  As  crystal  growth  progresses,  the 
crystal  solidifies  and  the  submerged  heater  is  slowly 
raised  in  order  to  maintain  a  constant  lower  melt  depth. 
Meyer  and  Ostrogorsky  (1997)  performed  experimental 
studies  of  the  flow  visualization  with  rotation  during 
the  Bridgman  process  with  a  submerged  heater.  The 
lower  melt  is  continuously  replenished  with  liquid  from 
the  upper  melt  with  a  composition  chosen  to  offset  the 
increasing  dopant  level  in  the  lower  melt  due  to  rejection 
at  the  crystal  growth  front  for  dopants  with  segregation 
coefficients  less  than  unity.  Unfortunately,  there  can  be 
severe  long-range  axial  segregation  in  the  crystal  if  the 
process  does  not  have  a  way  to  replenish  the  melt  with 
a  composition  that  offsets  rejection  of  species  along 
the  crystal-melt  interface  (Farrell  and  Ma,  2004).  Ostro¬ 
gorsky  and  Muller  (1994)  have  shown  that  this  modified 
bottom-seeded  VGF  method  produces  crystals  which 
exhibit  much  lower  defect  densities  and  much  more  dop¬ 
ant  uniformity  than  crystals  grown  by  the  LEC  process. 

Since  molten  gallium-antimonide  is  a  good  electrical 
conductor,  a  radial  electric  current  in  the  melt  can  be 
used  together  with  an  axial  magnetic  field  to  stir  the  melt 
in  order  to  control  the  dopant  distribution  in  the  semi¬ 
conductor  crystal,  which  depends  on  the  convective 
and  diffusive  transport  of  the  dopant  in  the  melt.  Wang 
et  al.  (2005)  have  demonstrated  that  a  radial  electric  cur¬ 
rent  of  a  few  amps  in  the  presence  of  a  weak  axial  mag¬ 
netic  field  can  produce  significant  stirring  rates  in  a 
gallium-antimonide  melt  having  a  50  mm  diameter  and 
1  cm  depth.  Some  of  the  advantages  of  electromagnetic 
(EM)  stirring  for  bottom-seeded  crystal  growth  are 
enhancement  of  the  radial  dopant  uniformity,  reduction 
of  thermal  stresses,  minimizing  of  defect  densities,  and 
the  ability  to  grow  at  faster  rates. 

The  large  number  of  adjustable  crystal  growth 
parameters  and  the  larger  number  of  ways  to  combine 
them  makes  process  optimization  through  trial-and- 
error  extremely  tedious.  Therefore,  models  that  predict 
melt  motion  and  heat  transfer  would  be  very  useful  for 
tailoring  the  melt  motion  in  order  to  optimize  the  com¬ 
plex  growth  process  in  a  timely  fashion.  In  the  present 
paper,  we  present  a  numerical  model  that  treats  the  melt 
motion  during  the  vertical  gradient  freeze  (VGF)  pro¬ 
cess  using  a  rotating  submerged  heater  and  with  a  uni¬ 
form,  steady  axial  magnetic  field  and  a  steady  radial 
electric  current. 


2.  Melt  motion 

This  paper  treats  the  melt  motion  of  selenium  in  a 
gallium-antimonide  melt  during  the  vertical  gradient 
freeze  crystal  growth  using  a  rotating  submerged  heater 
with  an  externally  applied,  uniform,  steady,  axial  mag¬ 
netic  field  Bz  combined  with  a  radial  electric  current. 
Here,  B  is  the  magnetic  flux  density,  while  r,  6  and  z 


are  the  unit  vectors  for  the  cylindrical  coordinate  sys¬ 
tem.  Our  dimensionless  problem  is  sketched  in  Fig.  1. 
The  coordinates  and  lengths  are  normalized  by  the  outer 
electrode’s  inner  radius  R ,  which  is  equal  to  the  crystal 
radius,  so  that  yQ  is  the  dimensionless  radius  of  the 
center  electrode,  y h  is  the  dimensionless  radius  of 
the  submerged  fused-silica  heater  housing,  and  b  is  the 
dimensionless  depth  of  the  lower  melt.  The  electrodes 
are  made  of  AXF-5Q  graphite,  which  is  manufactured 
by  Poco  Graphite.  This  graphite  has  an  electrical  con¬ 
ductivity  equal  to  1.47  x  1011  S/m  (Wang  et  al.,  2005). 
A  single  crystal  seed,  which  initiates  solidification,  lies 
at  the  bottom  of  the  fused  silica  crucible.  A  graphite  disc 
and  a  boron  nitride  disc  lie  below  the  crucible.  These 
discs  are  cooled  by  a  water-cooled  hearth  which  removes 
heat  along  the  bottom  of  the  crucible.  The  power  to  the 
submerged  heater  above  the  lower  melt  and  to  the  heat¬ 
ers  adjacent  to  the  periphery  of  the  crucible  are  adjusted 
so  that  the  crystal-melt  interface  is  nearly  planar,  and 
the  melt  adjacent  to  the  submerged  heater  is  held  at  a 
uniform  and  constant  temperature  Th.  As  the  crystal 
at  z  =  oct  solidifies,  the  crystal-melt  interface  moves  axi¬ 
ally  upward  at  the  dimensionless  rate  a  =  Ug/Uc ,  where 
Ug  is  the  constant  growth  rate  and  Uc  is  the  characteris¬ 
tic  velocity  of  the  melt  while  z  is  normalized  by  the  crys¬ 
tal’s  radius  R  and  t  is  time  normalized  by  RIUC.  The 
lower  melt  is  replenished  with  liquid  supplied  from  an 
upper  melt  to  a  small  annular  gap  between  the  sub¬ 
merged  heater  and  the  outer  electrode.  The  submerged 
heater  at  z  =  at  +  b  moves  upward  at  the  same  rate  so 


Center 


Outer 

Electrode 


Fused-Silica 

Crucible 
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r 


Bz 


Fig.  1.  VGF  process  using  submerged  heater  growth  with  a  uniform, 
steady,  axial  magnetic  field  Bz  and  with  a  steady  radial  electric  current 
where  coordinates  are  normalized  by  the  outer  electrode’s  inner  radius. 
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that  the  depth  of  the  melt  is  constant  throughout 
growth.  The  heater  is  rotated  at  angular  velocity  co, 
which  is  normalized  by  UJR.  A  long  cylindrical  graphite 
electrode  at  the  centerline  moves  with  the  submerged 
heater  and  carries  an  axially  downward  constant  DC 
electric  current.  The  electrical  conductivity  of  the  graph¬ 
ite  electrodes  is  much  larger  than  that  of  either  the  mol¬ 
ten  semiconductor,  the  fused- silica  heater  housing  or  the 
crystal,  so  that  the  electric  current  travels  through  the 
long  outer  graphite  electrode  tube  which  rests  inside 
the  periphery  of  the  crucible. 

A  magnetic  field  is  applied  with  a  copper  coil  outside 
the  heater  insulation,  while  an  electric  current  is  applied 
through  the  electrode  so  that  the  electromagnetic  (EM) 
body  force  created  by  this  electric  current  is  balanced 
by  the  inertial  force.  This  balance  gives  a  characteristic 
velocity  (Ma  et  al.,  2003), 


(i) 


where  R  is  the  inner  radius  of  the  outer  electrode  or  the 
radius  of  the  crystal,  and  p0  is  the  density  of  the  molten 
semiconductor  at  the  melting  temperature  Tm.  Here,  the 
characteristic  electric  current  density  is  Jc  =  II(2nr0bR ), 
where  r0  is  a  radial  position  that  can  be  used  to  provide 
an  approximate  estimate  of  the  radial  electric  current 
density.  For  a  crystal  grown  with  R  =  25  mm,  b  -  0.4 
and  r0  =  1  cm,  and  for  electromagnetic  stirring  produced 
with  i?  =  0.004  T  and  7=4  A,  our  estimation  of  Jc 
is  6366.2  A/m2  and  the  characteristic  velocity  is  Uc  = 
0.01027  m/s. 

The  electric  current  in  the  melt  produces  an  induced 
magnetic  field  which  is  superimposed  upon  the  applied 
magnetic  field  produced  by  the  external  magnet.  The 
characteristic  ratio  of  the  induced  to  applied  magnetic 
field  strengths  is  the  magnetic  Reynolds  number, 
Rm  -  HpcrUcR,  where  /ip  is  the  magnetic  permeability  of 
the  melt  and  o  is  the  electrical  conductivity  of  the  melt. 
For  all  crystal  growth  processes,  Rm  <C  1  and  the  addi¬ 
tional  magnetic  fields  produced  by  the  electric  currents 
in  the  melt  are  negligible. 

Ohm’s  law  is 


j  =  —V(j)+Nv  x  z,  (2) 

where  j  is  the  electric  current  density  normalized  by  /c,  </> 
is  the  electric  potential  normalized  by  JcRI(j,  and 
v  =  vri  +  ve0  +  vzz  is  the  melt  velocity  normalized  by 
Uc.  Here,  the  ratio  of  the  induced  electric  field  v  x  z  to 
the  static  electric  field  V</>  is  the  interaction  parameter, 
N  =  (jB2R/p0Uc .  For  the  present  process  with  B  -  0.004 
T  and  7=4  A,  N -  0.006  so  that  the  induced  electric 
field  is  negligible.  For  our  axisymmetric  problem,  there 
is  no  azimuthal  flow  of  electric  current,  and  the  electric 
potential  is  governed  by 

v2(p  =  0.  (3) 


We  treat  the  electrodes  as  perfect  conductors  which 
carry  a  uniform  electric  potential.  The  boundary  condi¬ 
tions  are 

(/)  =  0,  at  r  =  1,  for  —  1  ^  £  <  +1,  (4a) 


</>  =  at  C  =  +1,  for  0  ^r^ye.  (4b) 


Here,  £  =  —  1  +  2 (z  —  (xt)lb  is  a  rescaled  axial  coordi¬ 
nate  so  that  —  1  ^  £  <  +1.  The  fused-silica  heater,  the 
annular  gap  between  the  heater  and  the  outer  electrode, 
and  the  crystal-melt  interface  are  treated  as  electrical 
insulators  so  that  n  •  V</>  =  0,  where  n  is  the  outward 
unit  normal  vector.  The  electric  potential  is  given  by  a 
separation-of-variables  solution, 


<f>  =  y^„<Pfr/o(A„r) 

n= 0 


cosh(^(C+l)) 
cosh(2„Z?)  ’ 


(5) 


where  /0  is  the  Bessel  function  of  the  first  kind  and  zer¬ 
oth  order  and  Xn  are  the  eigenvalues  of  Jo-  Based  on  our 
choice  for  the  characteristic  electric  current  density  Jc , 
the  integral  of  the  radial  electric  current  over  the  melt 
depth  must  equal  b  at  r  -  r0,  so  that 


&o  —  b 


oo 

^  ^AnJ it(2nro) 

n= 1 


sinh(2„Z?) 

cosh(2„Z?) 


(6) 


where  7)  is  the  Bessel  function  of  the  first  kind  and  first 
order.  A  Galerkin  method  is  implemented  to  avoid  a 
Gibb’s  phenomena  associated  with  the  discontinuous 
boundary  conditions  along  £  =  +1.  The  electric  current 
originates  from  the  center  electrode  and  flows  radially 
outward  and  axially  downward,  as  reflected  in  the  con¬ 
tours  of  the  electric  potential  which  are  presented  in 
Fig.  2  for  ye  =  0.1905  and  b  -  0.4  for  which  &0  =  1.353. 

In  a  frame  of  reference  moving  with  the  fused-silica 
heater  and  crystal-melt  interface,  the  melt  motion  is 


Fig.  2.  Contours  of  the  electric  potential  in  the  melt  <p(r,  '()  for 
ye  =  0.1905  and  b  =  0.4. 
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steady.  With  the  Boussinesq  approximation,  the  equa¬ 
tions  governing  the  steady,  axisymmetric  melt  velocity 
and  temperature  are 

(▼•  V)v=  -S7p  +  ATz  +  ^-6  +  Re-'V2\,  (7a) 

V  •  v  =  0,  (7b) 

Pr  Re(\  •  V)r  =  V2r,  (7c) 

where  p  is  the  deviation  of  the  dimensional  pressure 
from  the  hydrostatic  pressure  for  a  uniform  density  nor¬ 
malized  by  p0t/2,  and  T  is  the  deviation  of  the  dimen¬ 
sional  temperature  from  the  melting  temperature  Tm 
normalized  by  ( Th  -  Tm).  Here,  the  Prandtl  number  is 
Pr  =  pcplk  while  the  Reynolds  number  is  Re  =  poUcRI 
p,  where  p,  cp  and  k  are  the  dynamic  viscosity,  specific 
heat  and  thermal  conductivity  of  the  melt,  respectively. 
The  Reynolds  number  is  related  to  the  strength  of  the 
electromagnetic  stirring,  i.e.,  Re  is  directly  proportional 
to  (BI)1/2.  In  Eq.  (7a),  the  ratio  of  the  buoyancy  force  to 
the  electromagnetic  body  force  is  A  -  pogf3(Th  —  Tm)/ 
JCB ,  which  we  refer  to  as  the  electromagnetic  (EM) 
buoyancy  parameter.  Here,  g  =  9.81  m2/s  and  /?  is  the 
thermal  volumetric  expansion  coefficient. 

We  use  a  frame  of  reference  moving  with  the  heater 
and  the  crystal-melt  interface,  so  the  no-slip  and  no¬ 
penetration  conditions  on  the  surface  of  the  outer 
electrode  are 


Vr 

=  o, 

at  r  =  1 , 

for  -  1  <  C  <  +1, 

(8a) 

ve 

=  0, 

at  r  =  1 , 

for  -  l  +1, 

(8b) 

vz 

=  —a, 

at  r  =  1 , 

for  —  1  <  £  <  +1. 

(8c) 

The  boundary  conditions  on  the  crystal-melt  inter¬ 
face  and  on  the  surfaces  of  the  center  electrode,  fused- 
silica  heater  and  the  annular  gap  between  the  heater 
and  the  crucible  where  the  lower  melt  is  replenished  are 


vr  =  0,  at  £  =  —  1,  for  0  ^  r  ^  1,  (9a) 

v0  =  0,  at  C  =  —  1 ,  for  0  <  r  ^  1,  (9b) 

vz  =  —a,  at  C  =  —  1 ,  for  0  <  r  <  1,  (9c) 

vr  =  0,  at  C  =  +1,  for  0  ^  r  ^  yh,  (9d) 

vq  =  cor,  at  C  =  +1,  for  o  <  r  <  7a,  (9e) 


Ve=J~y2  at^=+1’  f0r  7h  <  ^  1  > 

(9f) 

dz  =  0,  at  £  =  +1,  forO<r^yA,  (9g) 

vz=f{r),  at  C  =  +1,  for  yh  <  r  <  1.  (9h) 


We  use  simple  analytical  solutions  to  approximate 
the  azimuthal  and  axial  flows  of  the  replenished  fluid 
in  the  gap.  In.  Eq.  (9f),  we  approximate  the  azimuthal 
velocity  with  a  Taylor-Couette  flow  which  satisfies  the 
no-slip  conditions  given  by  Eqs.  (8b)  and  (9e).  In  Eq. 
(9h),  the  axial  velocity  of  the  fluid  that  is  replenished 
through  the  gap  is  approximated  with  a  Couette  flow. 
Assuming  that  the  density  of  the  crystal  and  melt  are 
the  same,  the  gap  between  the  heater  and  the  outer  elec¬ 
trode  replenishes  the  solidifying  melt  at  a  volumetric 
flowrate  noc,  so  that  /(r)  given  by 

fir)  =  «[-!+  Kf  1  -  r)  +  k2(  1  -  r2)\,  (10a) 


where 

_3(!+rD  „  _  2(1  +yh  +  y2h) 

(1-7/,)  (1  -7a)  (1  +7a) 


(10b,  10c) 


which  satisfies  the  no-slip  conditions  given  by  Eqs.  (8c) 
and  (9g). 

The  thermal  boundary  conditions  are 


T  =  0, 

aU  =  — 1, 

for  0  ^  r  ^  1, 

(1  la) 

T=  1, 

at  C  =  + 1 5 

for  0  ^  r  ^  1, 

(lib) 

3  T  „ 
37-°’ 

at  r  =  1 , 

for  -  1  <  C  ^  +1- 

(11c) 

We  have  assumed  that  the  rate  of  heat  transfer  from 
the  outer  electrode  to  the  melt  is  negligible. 

We  introduce  a  Stokes  stream  function  for  the  radial 
and  axial  velocities  in  the  meridional  circulations, 


l  0*A 

Vr 

r  oz 


1  # 

r  07’ 


(12a,  12b) 


which  identically  satisfies  conservation  of  mass  for  our 
axisymmetric  melt  motion.  For  given  values  of  A ,  Pr 
and  Re ,  we  solved  for  i //,  vq  and  T  using  a  Chebyshev 
spectral  collocation  method  with  Gauss-Lobatto  collo¬ 
cation  points  in  r  and  £.  Since  Eqs.  (7a)  and  (7c)  are 
non-linear,  we  used  a  Newton-Raphson  iterative  meth¬ 
od,  and  incrementally  increased  Re  from  zero.  After  the 
Newton-Raphson  method  converged  for  a  given  Re,  we 
used  a  continuation  method  to  obtain  initial  guesses  for 
the  coefficients  in  the  spectral  representations  for  (Re  + 
A  Re).  For  each  value  of  Re,  the  numbers  of  collocation 
points  was  varied  in  order  to  insure  that  the  results  were 
independent  of  this  numbers  and  to  insure  that  all 
boundary  layers  were  adequately  resolved.  For  example, 
for  Re  =  60  and  co  =  0.8542,  we  used  46  collocation 
points  in  the  radial  direction  and  41  collocation  points 
in  the  axial  direction.  When  we  doubled  the  collocation 
points  in  the  radial  direction,  the  minimum  value  of  the 
azimuthal  velocity  changed  by  0.03%  and  the  maximum 
value  of  azimuthal  velocity  changed  by  0.03%.  When  we 
doubled  the  collocation  points  in  the  axial  direction,  the 
minimum  and  maximum  values  of  the  azimuthal  veloc¬ 
ity  changed  by  only  0.0001%. 
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3.  Results 

We  investigate  the  effects  of  electromagnetic  stirring 
and  heater  rotation  on  the  melt  motion  with  yh  -  0.8, 
Pr  -  0.0442,  b  -  0.4,  and  Ug  =  2  mm/h  for  a  gallium- 
antimonide  melt.  The  parameters  as  a  function  of 
the  characteristic  velocity  are  Re  =  65,259.7  Uc,  Pet- 
2891.58  C/c,  A  =  0.0004697 1/“2,  and  a  =  5.556  x  10“7 
U~l,  where  Uc  is  in  m/s.  The  thermophysical  properties 
of  gallium-antimonide  are  presented  in  Table  1.  The 
characteristic  velocities  and  dimensionless  parameters 
are  summarized  in  Table  2. 

For  Re  -  60,  the  characteristic  velocity  is  Uc  = 
0.0009194  m/s.  The  dimensionless  parameters  are 
Pet  -  2.659,  A  -  555.7,  and  a  =  0.0006043.  We  present 
results  for  a  static  heater  with  co  =  0.  In  Fig.  3a-c,  we 
present  contours  of  the  azimuthal  velocity  vq(t,  Q,  the 
meridional  streamfunction  Q,  and  the  temperature 
T(r,  C),  respectively.  In  Fig.  3a,  the  electromagnetic 
(EM)  body  force  drives  a  clockwise  azimuthal  melt  mo¬ 
tion  and  the  minimum  and  maximum  values  of  vq  are 
—0.9248  and  0,  respectively.  This  minimum  value 
corresponds  to  a  dimensional  azimuthal  velocity  of 


Table  1 

Thermophysical  properties  for  molten  GaSb 


Property 

Value 

Density,  p0  (kg/m3) 

6030 

Dynamic  viscosity,  p  (Pas) 

0.00231 

Thermal  conductivity,  k  (W/mK) 

17.1 

Specific  heat,  cp  (J/kg  •  K) 

328 

Electrical  conductivity,  o  (S/m) 

1  x  106 

Thermal  coefficient  of  volumetric  expansion,  /?T  (K-1) 

0.0000958 

Table  2 

Characteristic  velocities  and  dimensionless  parameters 

Uc  =  0.0009194  m/s 

Uc  =  0.0004597  m/s 

Re 

60 

30 

Pet 

2.659 

1.3293 

A 

555.7 

2222 

a 

0.0006043 

0.001209 

-0.0008503  m/s.  In  Fig.  3b,  there  are  two  opposing  cir¬ 
culations  in  the  meridional  melt  motion  and  the  mini¬ 
mum  and  maximum  values  of  i j/(r,  £)  are  -0.002401 
and  0.001675,  respectively.  In  Fig.  3c,  we  can  see  that 


Fig.  3.  Contours  for  the  steady-state  melt  motion  with  Re  =  60  and  co  =  0:  (a)  azimuthal  velocity  ve(r,  '(),  (b)  meridional  streamfunction  i j/(r,  Q, 
(c)  temperature  T(r,  Q. 
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Fig.  4.  Contours  for  the  steady-state  melt  motion  with  Re  =  60  and 
co  =  —0.8542:  (a)  azimuthal  velocity  ve(r,  Q,  (b)  meridional  stream- 
function  Q. 

the  heat  transfer  is  dominated  by  conduction.  The  tem¬ 
perature  gradient  is  primarily  axial  with  a  slight  radial 
gradient. 

A  clockwise  rotation  of  the  heater  for  co  <  0  aug¬ 
ments  the  EM  stirring.  For  Re  -  60  and  co  =  —0.8542 
which  corresponds  to  a  rotation  rate  of  —0.3  rpm,  we 
present  contours  of  the  azimuthal  velocity  and  the 
meridional  streamfunction  in  Fig.  4a  and  b,  respectively. 
The  isotherms  are  very  similar  to  those  for  co  =  0  in  Fig. 
3c  in  which  the  heat  transfer  is  dominated  by  conduc¬ 
tion.  In  Fig.  4a,  the  heater’s  rotation  has  decreased  the 
minimum  value  of  ve  to  -1.0502  and  the  contours  are 
no  longer  closed.  This  minimum  value  corresponds  to 
a  dimensional  azimuthal  velocity  of  -0.0009656  m/s. 
In  Fig.  4b,  the  lower  circulation  takes  a  much  larger  vol¬ 
ume  of  the  melt  compared  with  the  co  =  0  case.  The  mag¬ 
nitude  of  the  radially  inward  flow  along  the  crystal-melt 
interface  has  increased.  The  minimum  and  maximum 
values  of  xj/  are  -0.003838  and  0.001125,  respectively. 


Fig.  5.  Contours  for  the  steady-state  melt  motion  with  Re  =  60  and 
co  —  —3.1417:  (a)  azimuthal  velocity  ve(r,  Q,  (b)  meridional  stream- 
function  Q. 

We  rotate  the  heater  more  quickly  at  co  =  -3.1417 
which  corresponds  to  a  rotation  rate  of  —  1 .2  rpm.  Con¬ 
tours  of  the  azimuthal  velocity  and  the  meridional 
streamfunction  are  presented  in  Fig.  5a  and  b,  respec¬ 
tively.  The  isotherms  are  very  similar  to  those  for 
co  =  0  in  Fig.  3c  in  which  the  heat  transfer  is  dominated 
by  conduction.  With  this  faster  rotation  rate,  the  mini¬ 
mum  value  of  vq  has  further  decreased  to  —2.694.  This 
minimum  value  corresponds  to  a  dimensional  azimuthal 
velocity  of  -0.002477  m/s.  The  increased  centrifugal 
body  force  is  reflected  in  Fig.  5a  where  the  contours  of 
the  azimuthal  velocity  are  nearly  dominated  by  the  flow 
due  to  rotation.  When  we  further  increased  the  rotation 
rate,  the  rotational  flow  dominated  the  EM  stirring. 
In  Fig.  5b,  the  minimum  and  maximum  values  of 
the  meridional  streamfunction  \j/  are  —0.01418  and 
0.0003021,  respectively.  The  lower  circulation  has  grown 
and  takes  the  entire  volume  of  the  melt,  in  which  the 
melt  flows  radially  inward  along  the  crystal-melt 
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Fig.  6.  Contours  for  the  steady-state  melt  motion  with  Re  =  60  and 
co  =  0.8542:  (a)  azimuthal  velocity  vo(r,  Q,  (b)  meridional  streamfunc- 
tion  0- 


Fig.  7.  Contours  for  the  steady-state  melt  motion  with  Re  =  60  and 
co  =  5.6950:  (a)  azimuthal  velocity  vo(r,  Q,  (b)  meridional  streamfunc- 
tion  Q. 


interface,  axially  upward  along  the  centerline,  radially 
outward  along  the  heater  and  axially  downward  along 
the  periphery  of  the  melt. 

For  Re  =  60,  we  investigate  a  counterclockwise  rota¬ 
tion  of  the  heater  for  which  co  >  0.  Because  the  EM  stir¬ 
ring  drives  a  clockwise  azimuthal  melt  motion,  a 
counterclockwise  rotation  of  the  heater  opposes  the 
EM  stirring.  For  co  =  0.8542,  Fig.  6a  and  b  show  the 
contours  of  the  azimuthal  velocity  and  the  meridional 
streamfunction,  respectively.  The  isotherms  are  very 
similar  to  those  for  co  =  0  in  Fig.  3c  where  the  heat  trans¬ 
fer  remains  dominated  by  conduction.  In  Fig.  6a,  the 
minimum  and  maximum  values  of  vq  are  —0.8120  and 
0.6735,  respectively.  These  minimum  and  maximum  val¬ 
ues  correspond  to  dimensional  values  of  -0.0007465  m/s 
and  0.0006192  m/s,  respectively.  In  Fig.  6b,  the  mini¬ 
mum  and  maximum  values  of  ij/  are  —0.001553  and 
0.001659,  respectively.  The  magnitude  of  the  radially  in¬ 


ward  flow  along  the  crystal-melt  interface  has  decreased 
compared  with  co  =  0. 

For  Re  -  60,  we  rotate  the  heater  more  quickly  at 
angular  velocity  co  =  5.6950  which  corresponds  to  a 
rotation  rate  of  2.0  rpm.  Contours  of  the  azimuthal 
velocity  and  the  meridional  streamfunction  are  pre¬ 
sented  in  Fig.  7a  and  b,  respectively.  The  isotherms 
are  very  similar  to  those  for  co  =  0  in  Fig.  3c  in  which 
the  heat  transfer  is  dominated  by  conduction.  With  the 
quicker  rotation  of  the  submerged  heater,  the  minimum 
and  maximum  values  of  vq  have  further  increased  to 
—0.4992  and  4.4903,  respectively.  These  minimum  and 
maximum  values  correspond  to  dimensional  values  of 
-0.0004590  m/s  and  0.004128  m/s,  respectively.  In  Fig. 
7b,  the  lower  circulation  takes  the  whole  volume  of 
the  melt.  The  minimum  and  maximum  values  of  x//  are 
—0.02381  and  0.0003021,  respectively.  The  minimum 
and  maximum  values  of  the  azimuthal  velocity  and 
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Table  3 


Minimum  and  maximum  values  of  vq  and  \j/ 


Re 

co*  (rpm) 

co 

^0,min 

vG,max 

*Amin 

*Amax 

60 

-1.2 

-3.1417 

-2.694 

0 

-0.01418 

0.0003021 

60 

-0.3 

-0.8542 

-1.0502 

0 

-0.003838 

0.001125 

60 

0 

0 

-0.9248 

0 

-0.002401 

0.001675 

60 

0.3 

0.8542 

-0.8120 

0.6735 

-0.001553 

0.001659 

60 

2.0 

5.6950 

-0.4992 

4.4903 

-0.02381 

0.0003021 

30 

-0.3 

-1.7085 

-1.3471 

0 

-0.001987 

0.0006043 

30 

0.3 

1.7085 

-0.2840 

1.3471 

-0.0009759 

0.0006043 

meridional  streamfunction  are  summarized  in  Table  3, 
where  co*  is  the  dimensional  rotation  rate  in  rpm. 

For  Re  -  30,  the  characteristic  velocity  is  Uc  = 
0.0004597  m/s,  and  the  dimensionless  parameters  are 
Pe,  =  1.3293,  A  -  2222,  and  a  =  0.001209.  For  a  clock¬ 
wise  rotation  co  =  —1.7085  for  which  the  rotation  rate 
is  -0.3  rpm,  we  present  the  contours  of  the  azimuthal 
velocity  and  meridional  streamfunction  in  Fig.  8a  and 


Fig.  8.  Contours  for  the  steady-state  melt  motion  with  Re  =  30  and 
co  —  —  1 .7085:  (a)  azimuthal  velocity  vd(r,  Q,  (b)  meridional  stream- 
function  Q. 


b  respectively.  The  isotherms  are  very  similar  to  those 
for  co  =  0  in  Fig.  3c  in  which  the  heat  transfer  is  domi¬ 
nated  by  conduction.  In  Fig.  8a,  the  minimum  and  max¬ 
imum  values  of  vq  are  —1.3471  and  0,  respectively.  This 
minimum  value  corresponds  to  a  dimensional  azimuthal 
velocity  of  —0.0006193  m/s.  The  clockwise  rotational 
flow  reinforces  the  clockwise  EM  stirring.  Com¬ 
pared  with  Re  =  60  and  a  rotation  rate  of  -0.3  rpm, 


Fig.  9.  Contours  for  the  steady-state  melt  motion  with  Re  =  30  and 
co  =  1.7085:  (a)  azimuthal  velocity  ve(r,  Q,  (b)  meridional  streamfunc¬ 
tion  Q. 
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the  contours  of  vq  in  Fig.  8a  reflect  a  much  stronger  rota¬ 
tional  flow  due  to  the  decreased  EM  stirring  for  Re  =  30. 
Because  the  magnitude  of  the  EM  stirring  is  propor¬ 
tional  to  (BI)1/2,  the  relative  importance  of  EM  stirring 
to  the  rotational  flow  has  decreased  significantly.  In 
Fig.  8b,  the  minimum  and  maximum  values  of  xj/  are 
-0.001987  and  0.0006043,  respectively. 

For  Re  -  30  with  a  counterclockwise  rotation  co  = 
1.7085  which  corresponds  to  a  rotation  rate  of  0.3 
rpm,  the  contours  of  the  azimuthal  velocity  and  meridi¬ 
onal  streamfunction  are  presented  in  Fig.  9a  and  b, 
respectively.  The  isotherms  are  very  similar  to  those 
for  co  =  0  in  Fig.  3c  in  which  the  heat  transfer  is  domi¬ 
nated  by  conduction.  In  Fig.  9a,  the  minimum  and  max¬ 
imum  values  of  vq  are  —0.2840  and  1.3471,  respectively. 
These  minimum  and  maximum  values  correspond  to 
dimensional  values  of  —0.0001306  m/s  and  0.0006192 
m/s,  respectively.  In  Fig.  9b,  the  minimum  and  maxi¬ 
mum  values  of  i//  are  -0.0009759  and  0.0006043,  respec¬ 
tively.  Compared  with  Re  =  60  and  a  rotation  rate  of  0.3 
rpm,  the  rotational  flow  more  strongly  opposes  the  EM 
stirring.  The  minimum  and  maximum  values  of  the  azi¬ 
muthal  velocity  and  meridional  streamfunction  are  sum¬ 
marized  in  Table  3. 

4.  Conclusions 

For  vertical  gradient  freezing  growth  using  a  sub¬ 
merged  heater  with  weak  magnetic  and  electric  fields, 
the  application  of  heater  rotation  can  have  a  strong 
influence  on  the  melt  motion.  A  larger  radially  inward 
flow  adjacent  to  the  crystal-melt  interface  is  desirable 
because  it  would  promote  compositional  uniformity  in 
the  crystal.  A  clockwise  heater  rotation  augments  the 
electromagnetic  (EM)  stirring  while  a  counterclockwise 
heater  rotation  opposes  the  EM  stirring.  The  magnitude 
of  the  radially  inward  flow  adjacent  to  the  crystal-melt 
interface  increases  as  the  Reynolds  number  (or  EM  stir¬ 
ring)  increases  and  as  the  rotation  of  the  heater  in  the 
clockwise  direction  increases. 
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